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Abstract 

In this lecture the next-to-leading order (NLO) corrections to the QCD Pomeron inter- 
cept obtained from the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation are discussed. It 
is shown that the BFKL Pomeron intercept when evaluated in non-Abelian physical renor- 
malization schemes with Brodsky-Lepage-Mackenzie (BLM) optimal scale setting does not 
exhibit the serious problems encountered in the modified minimal subtraction (MS) scheme. 
The results obtained provide an opportunity for applications of the NLO BFKL resumma- 
tion to high-energy phenomenology. One of such applications for virtual gamma-gamma 
total cross section shows a good agreement with preliminary data at CERN LEP. 
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1 Motivation 



The discovery of rapidly increasing structure functions in deep inelastic scattering (DIS) at 
HERA [|T| at small-x is in agreement with the expectations of the QCD high-energy limit. 
The Balitsky-Fadin-Kuraev-Lipatov (BFKL) 0, |^ resummation of energy logarithms is an- 
ticipated to be an important tool for exploring this limit. The leading order (LO) BFKL cal- 
culations predict a steep rise of QCD cross sections. Namely, the highest eigenvalue, cu'""^, 
of the BFKL equation |0] is related to the intercept of the Pomeron which in turn governs 
the high-energy asymptotics of the cross sections: a ~ s'^j'-i = gf^"""" The BFKL Pomeron 
intercept in the LO turns out to be rather large: ajp — 1 = uj^°-^ = 12 In 2 (a^/vr) ^ 0.55 for 
as = 0.2; hence, it is very important to know the next-to- leading order (NLO) corrections. 
In addition, the LO BFKL calculations have restricted phenomenological applications be- 
cause, e.g., the running of the QCD coupling constant as is not included and the kinematic 
range of validity of LO BFKL is not known. 

Recently the NLO corrections to the BFKL resummation of energy logarithms were 
calculated; see Refs. ^ and references therein. The NLO corrections |^ to the highest 
eigenvalue of the BFKL equation turn out to be negative and even larger than the LO 
contribution for as > 0.157. At such circumstances the phenomenological significance of the 
NLO BFKL calculations seems to be rather obscure. 

However, one should stress that the NLO calculations, as any finite-order perturbative 
results, contain both renormalization scheme and renormalization scale ambiguities. The 
NLO BFKL calculations 0, ^ were performed by employing the modified minimal subtrac- 
tion scheme (MS) [^] to regulate the ultraviolet divergences with arbitrary scale setting. 

In this work we consider the NLO BFKL resummation of energy logarithms 0, ^ in 
physical renormalization schemes in order to study the renormalization scheme dependence. 
To resolve the renormalization scale ambiguity we utilize Brodsky-Lepage-Mackenzie (BLM) 
optimal scale setting |^. We show that the reliability of QCD predictions for the intercept 
of the BFKL Pomeron at NLO when evaluated using BLM scale setting within non-Abelian 
physical schemes, such as the momentum space subtraction (MOM) scheme Q or the 
T-scheme based on T — > ggg decay, is significantly improved as compared to the MS- 
scheme. This provides a basis for applications of NLO BFKL resummation to high-energy 
phenomenology. Certain aspects of this work were presented in Ref. [ffO . 



2 BFKL in Physical Renormalization Schemes 

We begin with the representation of the MS- result of NLO BFKL 0, |^ in physical renormal- 
ization schemes. Although the MS-scheme is somewhat artificial and it lacks a clear physical 
picture, it can serve as a convenient intermediate renormalization scheme. The eigenvalue of 
the NLO BFKL equation at transferred momentum squared t = in the MS-scheme 0, 
can be represented as the action of the NLO BFKL kernel (averaged over azimuthal angle) 
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on the LO eigenfunctions {Q\IQ\ 
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= 2^(1)- ^(1/2 + ^(1/2 -zz/) 
is the function related with the LO eigenvalue, ip = F'/F denotes the Euler ip- function, the 



I/- variable is a conformal weight parameter |]TT|, Nc is the number of colors, and (5i.2 are 
the virtualities of the reggeized gluons. 

The calculations of Refs. [|, |^ allow us to decompose the NLO coefficient rj^ of Eq. 
(H) into /3-dependent and the conformal (/3-independent) parts: 
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Here /?o = (ll/3)A^'c — {2/3)Nf is the leading coefficient of the QCD /3-function, A^^ is the 
number of flavors, ({n) stands for the Riemann zeta-function, Li2(x) is the Euler dilogarithm 
(Spence-function). In Eq. (^) Np denotes flavor number of the Abelian part of the gg —>■ qq 
process contribution. The Abelian part is not associated with the running of the coupling 
1 12 1 and is consistent with the correspondent QED result for the 7*7* —>■ e^e~ cross section 



The /9-dependent NLO coefficient i^jjgi^) 



which is related to the running of the coupling, 
receives contributions from the gluon reggeization diagrams, from the virtual part of the one- 
gluon emission, from the real two-gluon emission, and from the non- Abelian part [|T2| of the 



gg — ^ qq process. There is an omitted term in r|^(z/) proportional to Xl^^) which originates 



from the asymmetric treatment of Qi and Q2, 
LO eigenfunctions [Q. 



it can be removed by the redefinition of the 
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The NLO BFKL Pomeron intercept then reads for A''( 



c 



a 



MS 



IP 

MS 



1=^M5W'>0) = 12 In 2 



^MsiQ"^ 



TT 



r— (0) 



20.12 - 0.1020iV^ + 0.06692/5, 
4 ^ -19.99 . 



1 + ^M5(0) 
) 



TT 



(6) 
(7) 



Physical renormahzation schemes provide small and physically meaningful perturbative 
coefficients by incorporating large corrections into the definition of the coupling constant. 
One of the most popular physical schemes is MOM-scheme ^, based on renormahzation 
of the triple-gluon vertex at some symmetric off-shell momentum. However, in the MOM- 
scheme the coupling constant is gauge- dependent already in the LO, and there are rather 
cumbersome technical difficulties. These difficulties can be avoided by performing calcula- 
tions in the intermediate MS-scheme, and then by making the transition to some physical 
scheme by a finite renormahzation [||. In order to eliminate the dependence on the gauge 
choice and other theoretical conventions, one can consider renormahzation schemes based 
on physical processes [0], e.g., V-scheme based on heavy quark potential. Alternatively, one 
can introduce a physical scheme based on the T ggg decay using the NLO calculations of 
Ref. ^ 
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A finite renormahzation due to the change of scheme can be accomplished by a transfor- 
mation of the QCD coupling M : 



as 



as 



1 + 

TT . 



where T is some function of Nc, Np and, for the MOM- scheme, of a gauge parameter ^. 
Then the NLO BFKL eigenvalue in the MOM-scheme can be represented as follows: 
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For the transition from the MS-scheme to the MOM-scheme the corresponding T-function 
has the following form 0]: 
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where I = -2 Jq dx \n{x)/[x^ - a; + 1] ~ 2.3439. 
Likewise, one can obtain for the V-scheme f? 
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Table 1: Scheme-transition function and the NLO BFKL coefficient in physical schemes 



Scheme 


Ji _ rpconf _j_ rpfi 


r(0) = r™"^(0) +r^(0) 


r(0) 
(iVf = 4) 


M 


e = o 


7.471 - 1.281/3o 


-12.64 - 0.1020A^F - 1.214/?o 


-22.76 





e = i 


8.247- 1.281/3o 


-11.87 - 0.1020iVi. - 1.214/?o 


-21.99 


M 




8.790 - 1.281/3o 


-11.33 - 0.1020A^F - 1.214/5o 


-21.44 


V 


2 - 0.4167/^0 


-18.12 - 0.1020A^F - 0.3497/5o 


-21.44 


T 


6.47 - 0.923/3o 


-13.6 - 0.102A^j. - 0.856/?o 


-21.7 



and, by the use of the results of Ref . [0] , for the T- scheme: 



Tt 



6.47 



2.77 



c 



(12) 



One can see from Table ^ that there is no a strong renormalization scheme dependence, 
though the problem of a large NLO BFKL coefficient remains. A large size of the perturbative 
corrections leads to a significant renormalization scale ambiguity. 



3 Optimal Renormalization Scale Setting 

The renormalization scale ambiguity problem can be resolved if one can optimize the choice 
of scales and renormalization schemes according to some sensible criteria. In the ELM opti- 
mal scale setting [j^, the renormalization scales are chosen such that all vacuum polarization 
effects from the QCD /3-function are resummed into the running couplings. The coefficients 
of the perturbative series are thus identical to the perturbative coefficients of the corre- 
sponding conformally invariant theory with (3 = 0. The BLM approach has an important 
advantage of resumming the large and strongly divergent terms in the perturbative QCD 
series which grow as nlfa^/^o]") ^-S-, the infrared renormalons associated with coupling con- 
stant renormalization. The renormalization scales in the BLM approach are physical in the 
sense that they reflect the mean virtuality of the gluon propagators 0. 

The BLM scale setting can be applied within any appropriate physical scheme. In the 
present case one can show that within the V-scheme (or the MS-scheme) the BLM procedure 
does not change significantly the value of the NLO coefficient r(z/). This can be understood 
since the V-scheme as well as MS-scheme are primarily adjusted to the case when, in the 
LO, there are dominant QED (Abelian) type contributions, whereas, in the BFKL case, the 
LO gluon-gluon (non- Abelian) interactions are important. 

Therefore, from the point of view of BLM scale setting, one can separate QCD processes 
into two classes specifying whether gluons are involved into the LO or not. Thus one can 
expect that in the BFKL case it is appropriate to use a physical scheme which is adjusted to 
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Table 2: The NLO- BFKL-Pomeron intercept in the BLM scale setting within non-Abelian 
physical schemes 



Scheme 


{Nf = 4) 




= 1 GeV 


= 15 GeV" 


= 100 GeV" 


M 

M 


e = o 


-13.05 


0.134 


0.155 


0.157 




-12.28 


0.152 


0.167 


0.166 


e = 3 


-11.74 


0.165 


0.175 


0.173 


T 


-14.01 


0.133 


0.146 


0.146 



non-Abelian interactions in the LO. One can choose the MOM- scheme based on the sym- 
metric triple-gluon vertex p, ^ or the T-scheme based on T — ggg decay. The importance 
of taking into account this circumstance for vacuum polarization effects can be seen from 
the "incorrect" sign of the /3o-term for rjjg in the unphysical MS-scheme (Eq. (|^)). 

Adopting BLM scale setting, the NLO BFKL eigenvalue in the MOM-scheme is 
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The /9-dependent part of the rMOM{^) defines the corresponding BLM optimal scale 
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Taking into account the fact that Xl{^) 



2^^ ^ 3 V 3 



-21n(z/) at z/ ^ oo, one obtains at large v 



Q¥Fm\^) = Q'-exp 



1 



(15) 



(16) 



At z/ = we have Qblm\0) = Q^(4exp[2(l + 2//3) - 5/3]) ~ QH27. Note that 
Qblm'^{^) contains a large factor, exp[— 4T^q^//3o] = exp[2(H-2//3)] ~ 168, which reflects 
a large kinematic difference between MOM- and MS- schemes [0, 0], even in an Abelian 
theory. 

Analogously, one can implement the BLM scale setting in the T-scheme (Table 0). 

Figures ^ and ^ give the results for the eigenvalue of the NLO BFKL kernel. We have 
used the QCD parameter A = 0.1 GeV which corresponds to as = 47r/[/5o ln((5^/A^)] ~ 0.2 
at = 15 GeV^. Also, the generalizations |^ [T^, |TB| of the /3 -function in the running 
coupling and of flavor number for continuous treatment of quark thresholds have been used. 

One can see from Fig. |l|, that the maximum which occurs at non-zero u is not as 
pronounced in the BLM approach compared to the MS-scheme, thus it does not serve as a 
good saddle point at high energies. 
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Figure 1: i/-dependence of the NLO BFKL eigenvalue at = 15 GeV^: BLM (in MOM- 
scheme) - solid, MOM-scheme (Yennie gauge: C = 3) - dashed, MS-scheme - dotted. LO 
BFKL {as = 0.2) - dash-dotted. 



One of the striking features of this analysis is that the NLO value for the intercept of 
the BFKL Pomeron, improved by the BLM procedure, has a very weak dependence on the 
gluon virtuality Q^. This agrees with the conventional Regge theory where one expects a 
universal intercept of the Pomeron without any Q^-dependence. The minor Q^-dependence 
obtained, on one side, provides near insensitivity of the results versus precise value of A, 
and, on the other side, leads to the approximate scale and conformal invariance. Thus one 
may use conformal symmetry |Tl|, ^ for the continuation of the present results to the 



case t 7^ 0. 

Therefore, by applying the BLM scale setting within non-Abelian physical schemes 



(MOM- and T- schemes), we do not face the serious problems [pT| , p^ , |23[| which were 
present in the MS-scheme, e.g., oscillatory cross section disbehavior based on the saddle 
point approximation , and a somewhat puzzling analytic structure |^ of the MS-scheme 
result i, |. 

Since the BFKL equation can be interpreted as a "quantization" of the renormalization 
group equation [|19] , it follows that the effective scale should depend on the BFKL eigenvalue 
uj, associated with the Lorentz spin, rather than on u. Thus, strictly speaking, one can use 
the above effective scales as function of u only in "quasi-classical" approximation at large Q^. 
However, the present remarkable Q^-stabilty leads us to expect that the results obtained with 
LO eigenfunctions may not change considerably for t ^ due to the approximate conformal 
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Figure 2: Q^-dependence of the BFKL Pomeron intercept in the NLO. The notation is as in 
Fie. 1. 



mvariance. 



4 Other Approaches to Perturbation 
Theory Optimization 

Now we consider briefly the NLO BFKL within other approaches to the optimization of per- 
turbative theory, namely, fast apparent convergence (FAG) ||2^ and the principle of minimal 
sensitivity (PMS) 

By the use of the FAC P^, one can obtain 



Q'faM = Q^exp 



as{QlAM) 



71 



Po 



(17) 
(18) 



In the MS-scheme at u = 0, ujfac = 0.33 - 0.26 for = 1 - 100 GeV^. However, the 
NLO coefficient r(z/) and hence FAC effective scale, each have a singularity at z/q — 0.6375 
due to a zero of the XL(z/)-function. 
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In the PMS approach p5[ the NLO BFKL eigenvalue reads as follows 



^^PA^s(Q^ ^) = NcxM 



apMs{Q'^{v)) 



TT 



{C /2)apMs/T^ 



(19) 



1 + CapMs/T^ 

where the PMS effective coupling apMs is ^ solution of the following transcendental equation 



^ , ry-, I CapMs/T^ 

+ C m' 



C/2 



OiPMS 



1 + CapMsl TT / 1 + CapMsl tt 



r(z/), 



(20) 



with C = /5i/(4/?o) and A 
upj^js = 0.23 - 0.20 for 



102 — 38Np/3. At u = one obtains in the MS-scheme 
1 — 100 GeV^ but, by the same reason as in the FAG case, 
the PMS effective coupling has a singularity at z/q. Thus, the application of the FAG and 
PMS scale setting approaches to the BFKL eigenvalue problem leads to difficulties with the 
conformal weight dependence, which is an essential ingredient of the BFKL calculations. 
The unphysical behavior of the FAG and PMS effective scales for jet production processes 



has been noted in Refs. 26 



The problem can be resolved by the expansion of xl('^) -function near its zero to avoid 
unphysical behavior of the optimization procedure. 



5 Application for Gamma- Gamma Scattering 



The gamma-gamma total cross section calculated with the resummation of the leading energy 
logarithms was considered in P, |2^ . 

The total cross section of two unpolarized gammas with virtualities Qa and Qp in the 



LO BFKL [g7|, I reads as follows: 



, u-rr T T^QaQb i 2.71 \ \QrJ I 



i,k=T,L 



(21) 



with the gamma impact factors in the LO for the transverse and longitudinal polarizations: 



IT 



q I 2 



r(2 - ii^)T{2 + iiy) 



(22) 



FM = acEnas [T.4y— r(2 - ..)r(2 + ..) ' ' ' ^^3) 

where Regge scale parameter sq is proportional to a hard scale ~ Q%Qb^ ^ being the 
Euler F-function, and is the quark electric charge. 



9 



In the NLO BFKL case one should obtain the formula analogous to LO BFKL (Eq. pT]) . 
It has been demonstrated in Ref. that the infrared singularities at the NLO are cancelled 
out for impact factors of colorless particles. Therefore, in the NLO both the kernel of the 
BFKL equation and impact factors are infrared safe which confirm a self-consistence of such 
factorization scheme. 



While exact NLO impact factors of gamma are not known yet ||30| one can use the LO 



impact factors of Eqs. (p2| - ^ ) |[T3| , ^, implying that the main NLO corrections come 



from the NLO BFKL subprocess rather than from the impact factors []3T| , |32[| . Thus, in 
the NLO BFKL one can have Eq. (^) but with uj{Q'^,v) taken in the NLO. To imply 
the BLM procedure to the total cross section one can see that one can imply the BLM 
procedure directly to the NLO BFKL eigenvalue uj{Q'^, v) within the accuracy up to the 
next-to-next-to-leading order (NNLO) and higher subleading terms. 




Figure 3: The NLO BFKL Pomeron vs preliminary L3 data on virtual gamma-gamma cross 
section (with subtracted quark-box contribution) at energy 91 GeV of the e^e~ collisions. 
Solid curves: NLO BFKL in BLM; dashed: LO BFKL, and dotted: LO contribution. Two 
different choices of the Regge scale: sq = and sq = 2Q^. 

For numerical calculations the NLO BFKL eigenvalue ci;(Q^, u) in the MOM-scheme 
(Yennie gauge: ^ = 3) has been used. 

In Figs. ^ - 1^ the comparison of BFKL predictions in the LO and NLO BFKL [^, ^ 



improved by the BLM procedure with L3 Collaboration data pBl EM from CERN LEP is 



shown. Different curves reflect uncertainty with the choice of the Regge scale parameter 
which indicates when the asymptotic regime starts. At infinite collision energies, the cross 
sections do not depend on this scale parameter sq. For present calculations, two variants have 
been choosen sq = and Sq = 2Q^, where for symmetric virtuality case = Q\ = Q%. 
One can see from Figs. |5 - |^ that the LO BFKL predictions overestimate the L3 data, 
while the agreement of the NLO BFKL improved by the BLM procedure is reasonably well, 
especially at higher energies of LEP2 ^Se+e- = 183 - 189 GeV (Figs. ^, |^). One can notice 
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also that sensitivity of the NLO BFKL results with respect to the Regge parameter sq is 
much smaller than in the case of the LO BFKL. Recent OPAL Collaboration data are 
also in a good agreement with the NLO BFKL predictions. 



L3 data from LEP2 
1 VS„^ = 1 83 GeV 
<Q'> = 14 GeV' 



NLO BFKL+BLM 

LO BFKL 

LO 




Y=log(S„/<Q'>) 



Figure 4: The same as Fig. 3, but for energy 183 GeV of e^e collisions. 




Figure 5: The same as Figs. 3,4, but for energy 189 GeV of e'^e collisions. 

The gamma-gamma scattering is attractive from viewpoint that it is theoretically more 
controllable rather than hadron-hadron and lepton-hadron collisions where non-perturbative 
hadronic structure functions are involved. In addition, in the gamma-gamma scattering 
the unitarization (screening) corrections due to multiple Pomeron exchange would be less 



important than in hadron collisions. It was shown in Refs. |3g] that the unitarization 



11 



corrections in hadron collisions can lead to higher value of the (bare) Pomeron intercept than 
the effective intercept value. Since the hadronic data fit yields about 1.1 for the effective 
intercept value |3^, then the bare Pomeron intercept value should be above this value. 
So that, in case of small unitarization corrections in the gamma- gamma scattering at large 
one can accomodate the NLO BFKL Pomeron intercept value 1.13-1.18 along with larger 
unitarization corrections in hadronic scattering where it can lead to a smaller effective 
Pomeron intercept value about 1.1 for hadronic collisions. 



6 Summary 



There have been a number of recent papers which analyze the NLO BFKL predictions in 
terms of rapidity correlations t-channel unitarity angle-ordering |Q, double 

transverse momentum logarithms p3| , |45| , 

m 



I and ELM scale setting for deep inelas- 
This requires a further study to find relations between such 



tic structure functions 
approaches. 

To summarize, we have shown that the NLO corrections to the BFKL equation for the 
QCD Pomeron become controllable and meaningful provided one uses physical renormaliza- 
tion schemes relevant to non-Abelian gauge theory. BLM optimal scale setting automat- 
ically sets the appropriate physical renormalization scale by absorbing the non-conformal 
/3-dependent coefficients. The strong renormalization scale dependence of the NLO correc- 
tions to BFKL resummation then largely disappears. This is in contrast to the unstable 
NLO results obtained in the conventional MS- scheme with arbitrary choice of renormaliza- 
tion scale. A striking feature of the NLO BFKL Pomeron intercept in the BLM approach is 
its very weak Q^-dependence, which provides approximate conformal invariance. The NLO 
BFKL application to the total gamma-gamma cross section shows a good agreement with 
the preliminary L3 data at the CERN LEP2 energies. The results presented here open new 
windows for applications of NLO BFKL resummation to the high-energy phenomenology. 
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Note added at Proof. After presentation of this lecture, the OPAL and L3 Collaborations 
at the CERN LEP accumulated statistics at higher energies and finalized their data 
The final OPAL and L3 data ||4^, although presented in a different way, show even better 
agreement with our earlier predictions [^l[ 
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